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FOREWORD 


The  finite  element  code  STAGS  is  used  to  obtain  buckling  load  predictions 
for  a  number  of  ring  reinforced  circular  cylindrical  shells  for  which  experi¬ 
mental  results  are  reported  in  the  open  literature.  However,  as  the  thickness 
of  pressure  hulls  increases  and  the  ratio  of  radius  to  thickness  decreases, 
failure  may  occur  by  plastic  yielding  and  axisymmetric  collapse  at  a  limit  point 
For  these  cases,  predictions  using  STAGS  are  substantially  less  successful.  The 
possible  existence  of  locked-in  residual  stresses  or  strains  is  another  source 
of  discrepancy.  Agreement  is  reasonable,  especially  in  the  elastic  range,  in 
view  of  the  difficulties  of  this  type  of  analysis. 

This  work  represents  three  work-months  of  effort  under  funding  of  the  Naval 
Sea  Systems  Command,  whose  support  is  gratefully  acknowledged. 
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Figure 


ILLUSTRATIONS 


FULL  CYLINDER  CIRCUMFERENTIALLY  AND  HALF  AXIALLY  DISPLAYING 
CIRCULAR  ENDS  4  AND  3  AND  ASSUMED  LONGITUDINAL  ENDS  2 

AND  4  AT  8-0  and  0-ir  RESPECTIVELY . 

MODE  1  FOR  MODEL  I,  HALF  MODEL  OF  BR-5  . 


MODE  I  FOR  MODEL  2,  FULL  MODEL  OF  BR-5 


MODE  1  FOR  MODEL  1A  (BR-5  INCLUDING  1  BAY  WITH  BULKHEAD) 
SUBJECT  TO  CONTINUITY  CONDITIONS  AND  W=0  AT  BULKHEAD  END. 

EXTERNAL  PRESSURE-DEFLECTION/THICKNESS  CURVE  OF  MODEL  3 
UP  TO  APPROXIMATE  STATIC  COLLAPSE  . 


MODE  1  FOR  MODEL  4  (ALUMINUM  MODEL  G  OF  MIDGLEY  AND  JOHNSON) 
SUBJECT  TO  SIMPLY  SUPPORTED  CONDITIONS  AT  END  I  . 


MODE  1  FOR  MODEL  5  (ALUMINUM  MODEL  G  OF  MIDGLEY  AND  JOHNSON) 
SUBJECT  TO  CONTINUITY  END  CONDITIONS  AND  W#0  AT  END  I  .  .  . 

MODE  1  FOR  MODEL  6  (ALUMINUM  MODEL  G  OF  MIDGLEY  AND  JOHNSON) 
SUBJECT  TO  CONTINUITY  CONDITIONS  AND  W=0  AT  END  1  . 

MODE  1  FOR  MODEL  8  (ALUMINUM  MODEL  F  OF  MIDGLEY  AND  JOHNSON) 
SUBJECT  TO  SIMPLY  SUPPORTED  CONDITIONS  AT  END  I  . 


MODE  1  FOR  MODEL  9  (ALUMINUM  MODEL  F  OF  MIDGLEY  AND  JOHNSON) 
SUBJECT  TO  CONTINUITY  CONDITIONS  AND  W^O  WITHOUT  END 
RING  AT  END  1  . 

MODE  I  FOR  MODEL  10  (ALUMINUM  MODEL  F  OF  MIDGLEY  AND  JOHNSON) 
SUBJECT  TO  CONTINUITY  CONDITIONS  AND  W^O,  WITH  END  RING 
AT  END  I  . 

MODE  1  FOR  MODEL  II  (ALUMINUM  MODEL  F  OF  MIDGLEY  AND  JOHNSON) 
SUBJECT  TO  CONTINUITY  CONDITIONS  AND  W=0  AT  END  1  . 

MODE  1  FOR  MODEL  12  (STEEL  1/5  SCALE  MODEL  OF  KINRA)  SUBJECT 
TO  CONTINUITY  CONDITIONS  AND  W-0  AT  END  1  . 
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Figure 


ILLUSTRATIONS  ( Cent . ; 


MODE  1  FOR  MODEL  13  (BR-4,  EXTERNALLY  STIFFENED  STEEL  MODEL/ 
SUBJECT  TO  CONTINUITY  CONDITIONS  AND  W=0  AT  END  1  ...  . 


DEFORMED  SHAPE  OF  MODEL  14  (STEEL  MODEL  BR-4  j  AT  500  LB/  IN*- 

MODE  1  FOR  MODEL  17  (1.965  LONGER  THAN  MODEL  8)  SUBJECT  TO 
SIMPLY  SUPPORTED  CONDITIONS  AT  END  I  . 

MODE  1  FOR  MODEL  18  (1.965  LONGER  THAN  MODEL  9)  SUBJECT  TO 
CONTINUITY  CONDITIONS  AND  W#0,  WITHOUT  END 
RING  AT  END  I  . 

MODE  1  FOR  MODEL  19  (1.965  LONGER  THAN  MODEL  10)  SUBJECT  TO 
CONTINUITY  CONDITIONS  AND  W^O,  WITH  END  RING  AT  END  1  .  . 

MODE  1  FOR  MODEL  20  (1.965  LONGER  THAN  MODEL  11)  SUBJECT  TO 
CONTINUITY  CONDITIONS  AND  W“0  AT  END  1  . 
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TEST  MODEL  STRUCTURAL  DETAILS  . 

TEST  MODEL  MATERIAL  DETAILS  . 

ACTUAL  DIMENSIONS  USED  IN  COMPUTATIONS  . 

MODELING  DETAILS  . 

COMPARISON  BETWEEN  FINITE  ELEMENT  ANALYSIS  PREDICTIONS 
AND  EXPERIMENTAL  RESULTS  . 
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CHAPTER  1 


INTRODUCTION 


The  ultimate  objective  of  this  work  is  to  establish  reliable  finite  element 
procedures  for  predicting  the  effect  of  structural  damage,  for  example  due  to  a 
shock  wave,  on  the  collapse  strength  of  submarine  hulls.  However,  it  is  first 
necessary  to  demonstrate  the  validity  of  finite  element  calculations  for 
undamaged  structures.  Unfortunately,  the  open  literature  provides  little  infor¬ 
mation  on  finite  element  code  validation  for  buckling  of  stiffened  structures. 

l  ?  •  f* 

In  two  previous  studies  by  the  current  author,1*^  predictions  using  the 
code  STAGS'1  were  compared  to  those  from  other  methods  in  the  open  literature. 

The  present  report  compares  STAGS  predictions  to  published  experimental  results 
from  seven  different  tests.  The  agreement  is  reasonable,  particularly  in  view 
of  the  difficulties  of  this  type  of  analysis. 


^-Moussouros ,  M. ,  Comparisons  of  Static  Collapse  Pressure  Predictions  of  a 
Ring-Stiffened  Cylindrical  Shell  Subject  to  Hydrostatic  Pressure,  NSWC  TR 

81- 325,  3  Mar  1982. 

^Moussouros,  M. ,  Further  Results  on  the  Predictions  of  Collapse  Pressure  of 
a  Ring-Stiffened  Cylindrical  Shell  Subject  to  Hydrostatic  Pressure,  NSWC  TR 

82- 172,  Sep  1982. 

^Almroth,  B.  0.,  Brogan,  F.  A.,  and  Stanley,  G.  M. ,  Structural  Analysis  of 
General  Shells,  Vol.  II,  User  Instruction  for  STAGSC,  LMSC-D633873,  Apr  1979. 


1 


CHAPTER  2 


ANALYSIS 


A  number  of  papers  in  the  open  literature  contain  experimental  collapse 
pressure  values  for  ring-stif fenen  cylindrical  shells. 


^Slankard,  R.  C. ,  and  Nash,  W.  A.,  Tests  of  the  Elastic  Stability  of  a  Ring- 
Stiffened  Cylindrical  Shell,  Model  BR-5  U *1.705,  Subjected  to  Hydrostatic 
Pressure,  DTMB  Report  822,  May  1953. 

5siankard,  R.  C.,  Tests  of  Elastic  Stability  of  a  Ring-Stiffened  Cylindrical 
Shell,  Mooel  BR-4  (X»1.103)  Subjected  to  Hydrostatic  Pressure,  DTMB  Report 
876,  Feb  1955. 

^Kirstein,  A.  F. ,  and  Slankard,  R.  C. ,  An  Experimental  Investigation  of  the 
Shell-Instability  Strength  of  a  Machined,  Ring-Stiffened  Cylindrical  Shell 
Under  Hydrostatic  Pressure  (Model  BR-4A) ,  DTMB  Report  997,  Apr  1956. 

'Lunchick,  M. ,  and  Overby,  J.  A.,  An  Experimental  Investigation  of  the  Yield 
Strength  of  a  Machined  Ring-Stiffened  Cylindrical  Shell  (Model  BR-7M)  Under 
Hydrostatic  Pressure,  DTMB  Report  1255,  Nov  1958. 

^DeHart,  R. ,  and  Basdekas ,  N.  L. ,  "Investigation  of  Yield  Collapse  of 
Stiffened  Circular  Cylindrical  Shells  with  a  Given  Ou t -of -Roundness , "  in 
Collected  Papers  on  Instability  of  Shell  Structures-1962,  NASA  TN  D-1510, 
1962,  pp.  245-253. 

^Midgley,  W.  R. ,  and  Johnson,  A.  E.,  Jr.,  "Experimental  Buckling  of  Internal 
Integral  Ring-Stiffened  Cylinders,"  Experimental  Mechanics,  Jul  1973,  pp - 
145-153. 

^Kinia,  R.  K.,  "Hydrostatic  and  Axial  Collapse  Tests  of  Stiffened  Cylinders," 
Paper  2685  in  Offshore  Technology  Conference,  1976,  pp.  765-788. 

^Galletly,  G.  D. ,  Slankard,  R.  C.,  and  Wenk,  E. ,  Jr.,  "General  Instability 
of  Ring-Stiffened  Cylindrical  Shells  Subject  to  External  Hydrostatic 
Pressure — A  Comparison  of  Theory,"  Journal  of  Applied  Mechanics,  Vol.  25, 

No.  2,  Jun  1958,  pp.  259-266. 

^Reynolds.  T.  E.,  and  Blumenberg,  W.  F.,  General  Instability  of  Ring-Stif fenea 
Shells  Subject  to  External  Hydrostatic  Pressure,  DTMB  Report  1324,  Jun  1939. 
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Several  of  these  studies  ‘  provide  enough  data  to  permit  moueiing 
(e.g.,  mo<^el  design,  material  properties).  Unfortunately,  in  our  judgment,  the 
remainder*-4-^9  do  not  give  sufficient  information,  out  they  are  listen  here 
for  completeness. 

Table  1  gives  the  structural  details  of  the  models  to  be  examined,  wnile 
Table  2  gives  material  properties. 

The  various  models  have  been  analyzed  using  the  general  purpose  finice 
element  code  STAGS.  Tables  3  and  4  show  details  of  the  finite  element  analysis 
including  the  dimensions  actually  used  in  the  computat ions ,  the  number  of 


l^Krenzke,  m.  A.,  Effect  of  Initial  Deflections  ana  Residual  Welding  Stresses 
of  Elastic  Behavior  and  Collapse  Pressure  of  Stiffened  Cylinders  Subjected  to 
External  Hydrostatic  Pressure,  DTMB  Report  1327,  Apr  1960. 

l^Blumenberg,  W.  F. ,  and  Reynolds,  T.  E.,  Elastic  Instability  of  Ring- 

Stiffened  Cylinders  with  Intermediate  Heavy  Frames  Under  External  Hydrostatic 
Pressure ,  DTMB  Report  1588,  Dec  1961. 

l^Blumenberg,  W.  F. ,  Hydrostatic  Pressure  Tests  to  Determine  the  Effect  of 
Varying  Degrees  of  End  Fixity  on  the  Elastic  General  Instability  Strength  of 
Ring-Stiffened  Cylindrical  Shells,  DTMB  Report  2361,  May  1967. 

l^Batista,  r.  c. ,  and  Croll,  J.  G.  A.,  "Simple  Buckling  for  Pressurized 

Cylinders,"  Journal  of  Engineering  Mechanics  EM5,  ASCE,  Oct  1982,  pp. 927-944. 

17Slankard,  DTMB  Report  822. 

l^Slankard,  DTMB  Report  876. 

l^Kirstein,  DTMB  Report  997. 

^Lunchick,  DTMB  Report  1255. 

21 DeHart,  NASA  TN  D-1510. 

22Midgley,  pp.  145-153. 

23Kinra,  pp.  765-788. 

2^Galletly ,  pp.  259-266. 

25Reynolds,  DTMB  Report  1324. 

2^Krenzke,  DTMB  Report  1327. 

27Blumenberg ,  DTMB  Report  1588. 

^®Blumenberg ,  DTMB  Report  2361. 

29Bastista,  pp.  927-944. 
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degrees  of  freedom  (D.O.F.),  Che  discribucion  of  noaes ,  and  Che  boundary 
conditions.  The  models  are  assumed  Co  be  perfecClv  circular.  They  are 
discretized  using  fiac  place  quadri  lacera  1  elemenCs,  specifically  cne  STAGS  UJ 
elemenC.  Rigid  body  modes  are  removed  wich  appropriace  consCrainCs.  Axial 
loads  are  imposed  as  equivalenC  line  loads  along  Che  end  place  perimeCers.  This 
meChoa  neglecCs  local  bending  aC  end  place  inCerseccions .  ^Finally  Che  "aead 
pressure"  is  used.  The  base  load  pressure  used  in  1  lb/in“  and  Che  linear 
option  (linear  sCress  state)  of  STAGS  is  exercised,  unless  oCherwise.  raenCioned . 
Figure  1  illustrates  the  coordinate  system  and  the  notation  to  be  used  to  denote 
Che  displacements. 
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CHAPTER  3 

DISCUSSION  OF  RESULTS 


Table  5  displays  the  numerical  results  obtained  using  STAGS,  i.e.,  critical 
pressures  for  static  collapse,  buckling  mode,  type  of  boundary  conditions  used 
in  the  analysis,  and  type  of  stiffening.  Also  shown  are  experimental  collapse 
values  and  collapse  modes  where  available. 

Note  that  Models  1  and  2  give  a  buckling  pressure  of  106.3  lb/in-  for 
Model  1  and  102.6  lb/in^  for  Model  2  with  correspnding  oval  modes  (Figures  2 
and  3),  indicating  stiffener  collapse.*  Also  note  that  the  half  and  full  models 
circumferentially  do  not  give  identical  critical  pressures,  although  both  compu¬ 
tations  are  carried  out  using  the  same  base  load.  Model  1A,  which  incorporates 
an  additional  shorter  frame  spacing  with  a  rigid  bulkhead,  results  in  a  higher 
critical  pressure  of  155.596  lb/in^  and  a  lobar  mode  of  collapse  (Figure  41. 

For  Model  1A,  the  maximum  axial  and  hoop  stresses  are  -28.785  ksi  and  -35.631 
ksi,  respectively.  The  lower  experimental  pressure  (80  lb/in^)  is  attnbutea 
to  residual  strains  (model  was  not  stress  relieved)  rather  than  to  imperfections 
considering  that  in  the  case  of  ring-stiffened  cylinders,  the  rings  in  a  sense 
constitute  imperfections  far  larger  in  magnitude  than  the  ones  due  to  fabri¬ 
cation.  There  is  an  undeniaole  influence  due  to  imperfections,  however,  wh'ch 
according  to  CoppaJU  is  smoothed  out  with  increasing  pressure,  except  for  these 
imperfections  due  to  the  rings  themselves.  The  case  of  axial  compression  only 
is  excluded  in  this  discussion. 


A  nonlinear  incremental  analysis  was  performed  for  Model  3  (Figure  5).  It 
was  discontinued  at  a  pressure  9500  lb/in^  due  to  a  large  step  size  increment. 
This  is  comparable  to  the  experimental  collapse  pressure  of  9750  lbs/in'1'.  Note 
that  the  slope  of  the  pressure-aisplacement  curve  on  Figure  5  is  decreasing 
rapidly  in  this  vicinity. 


Model  4  (simply  supported  at  ena  1)  buckled  in  a  lol nr  pattern  (Figure  o) 
at  23.704  lb/in^,  while  Model  5  (Figure  7)  with  continuity  conditions  at  ena  1 
ovalized  at  21.262  lb/in^.  Model  6  with  W  =  0  at  end  1  (effect  of  bulkhead) 
buckled  at  23.669  lb/in^  in  a  diamond  shape  mode  between  rings  (Figure  8;. 


*In  Figure  2,  the  stiffener  at  midlength  is  shown  displaced  radially  inwards. 

-^Coppa,  A.  P. ,  "Measurement  of  Initial  Geometrical  Imperfections  of 

Cylindrical  Shells,"  AIAA  Journal,  Vol.  4,  No.  1,  Jan  1966,  pp.  172-175. 
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According  Co  Midgley  and  Johnson,  1  Che  stiffeners  were  overdesigned  co  iocalize 
the  instability  between  successive  rings.  Model  7  was  treated  using  an  incre¬ 
mental  nonlinear  analysis.  When  discontinued  at  31  lb/in~,  the  model  hau  not 
attained  the  collapse  pressure.  The  experimental  collapse  pattern  given  is  thac 
of  Midgley  and  Johnson. 32  Overall,  if  we  are  to  take  the  first  value  given 
(25  lb/in2),  the  finite  element  program  STAGS  gives  a  fairly  good  estimate  tor 
this  model.  The  experimental  model  represented  by  Models  4  through  7  was  a 
machined  structure. 

Models  9  and  10  (Figures  10  and  11)  yield  collapse  pressures  well  below  tne 
lowest  experimental  value  (65  in/in2)  reported.  Models  8  and  11  (Figures  9 
and  12)  failed  by  overall  instability  at  62.651  lb/in2  and  62.371  lb/in2, 
respectively.  These  models  were  based  on  Model  F.  J  The  ring  stiffeners  were 
underdesigned,  unable  to  prevent  overall  collapse  or  restrict  buckling  between 
them.  Figures  9  through  12  exhibited  overall  instability,  confirmed  by  Midgley 
and  Johnson.'*  An  oval  critical  mode  was  obtained  when  continuity  conditions 
at  end  1,  without  an  end  ring  (Model  9  of  Figure  10),  were  used  and  also  when  an 
end  ring  was  employed  (Model  10  of  Figure  11).  Since  the  stiffeners  were  not 
strong  enough,  they  collapsed  according  to  the  ring  formula,  which  in  this 
case  gives  5.266  lb/in2  (Table  5).  The  experimental  model  represented  by 
Models  8  through  11  was  a  machined  structure. 

On  the  other  hand,  Southwell's  formula,3**  treating  the  shell  as  fully 
unstiffened,  yields  about  6.60  in/in2  (this  is  not  included  in  Table  5)  as  the 
collapse  pressure.  Using  Brush  and  Almroth's  Figure  5.17  or  5.113^  with  L  = 
18.85  in.,  a  =  7.955  in.,  h  =  0.040  in.,  E  =  10.10  x  10**  lb/in2,  v  =  0.300, 

D  =  59.194  lb-in.,  Z  =  1065.8,  we  obtain  p~36  and,  therefore,  the  critical 
collapse  pressure  is  7.43  lb/in2.  It  is  evident  from  this  analysis  that  only 
a  very  long  cylinder  with  weak  stiffeners  should  buckle  as  displayed  on  Figures 
10  and  11.  On  the  other  hand,  when  the  stiffeners  are  strong  (Figure  7  of  Model 
5),  we  may  still  obtain  a  relatively  accurate  estimate  of  the  critical  pressure. 
It  is,  however,  associated  with  the  wrong  mode.  At  first  sight,  it  appears  that 
the  computed  collapse  pressures  for  Models  2  and  4  are  reasonably  good,  although 
the  predicted  collapse  modes  are  not.  It  is  considered  that  their  apparent  good 
agreement  should  not  be  trusted  too  much.  As  discussed  in  the  next  paragraph, 
the  lower  experimental  pressure  is  probably  due  to  residual  strains  and  initial 
imperfections  in  this  welded  cylindrical  model. 


^Midgley,  pp.  145-153. 

33Midgley,  pp.  145-153. 

33Midgley,  pp.  145-153. 

^Midgley,  pp.  145-153. 

33Moussouros ,  NSWC  TR  82-172. 

36Moussouros,  NSWC  TR  81-325. 

3^Brush,  D.  0.,  and  Almroth,  B.  0.,  Buckling  of  Bars,  Plates,  and  Shells, 
McGraw  Hill,  Inc.,  1975,  pp.  161-167. 


Model  12  (.Figure  clearly  shows  lobar  buckling  between  rings  a:  a  nigher 

pressure  (161.19  lb/in’~)  reported  on  Table  5  (110-115  lb/in-).  Another  unre¬ 
ported  model,  with  frame  spacing  twice  that  of  Model  12,  collapsed  at  78. oo 8 
lb/in2  in  the  lobar  mode.  The  experimental  moael  represented  by  Mooei  l-  was 
a  welded  structure. 


Model  13  (Figure  14)  collapsed  at  536.058  lb/in-  in  a  lobar  moae.  Ac 
this  pressure,  a  hoop  stress  of  -58.96  ksi  (higher  than  yield)  developed  ano  the 
linear . stress  state  must  be  abandoned.  Model  14  (Figure  15),  which ^exhibited 
plastic  deformation,  had  not  failed  up  to  a  pressure  of  512.5  lb/in^.  At  this 
point,  it  appears  that  residual  stresses  were  the  cause  of  this  large  discrepancy 
(390  in/in2).  A  model  identical  to  BR-4,  but  stress  relieved,  failed  at  550 
lb/in2  (compared  to  390  in/in2,  which  includes  residual  stresses). This  is 
close  to  the  reported  collapse  pressure  by  Bushnell^  and  the  BASOR  program 
(460  lb/in2). 


In  Model  15,  collapse  pressure  predictions  ignoring  plasticity  are  much 
higher  than  the  experimental  values  used  and  should  be  disregarded.  Model  16  is 
treated  by  nonlinear  analysis;  it  did  not  collapse  up  to  1500  lb/in2.  It  must 
be  stressed  at  this  point  that  except  for  Model  3,  for  which  the  stress  strain 
curve  was  given,  the  stress  strain  curve  had  to  be  approximated  and  adjusted, 
causing  another  source  of  differences. 


Models  17,  18,  19,  and  20  (Figures  16,  17,  18  and  19)  collapsed  in  a 
fashion  similar  to  Models  8,  9,  10,  and  11  respectively,  but  at  lower  critical 
pressures,  since  their  length  was  larger.  Observe  that  the  critical  collapse 
pressures  for  Models  18  and  19  (4.13  lb/in2  and  4.268  lb/in2)  using  continuity 
conditions  at  end  1  were  only  about  50%  of  the  critical  pressure  for  Models  9 
and  10  (7.257  lb/in2  and  7.524  lb/in2).  This  ratio  (0.569  to  0.567)  in  critical 
pressure  was  not  the  same  for  Models  17  anc;  20  (19.893  lb/in2  and  20.277 
lb/in2)  as  compared  to  Models  8  and  11  (62.651  lb/in2  and  62.371  lb/in2).  This 
last  range  of  ratios  in  critical  pressure,  as  4  is  reduced  to  0.50894,  becomes 
0.317  to  0.325,  suggesting  a  trend  similar  to  a  strut.  These  comments  follow, 
since  the  mode  here  is  that  of  general  instability  in  the  elastic  range.  It  is 
known  that,  when  the  length  of  a  strut  is  increased  from  4^  to  42,  the  critical 
pressure  within  the  elastic  range  is  reduced  to  the  ratio  (4^/42)  .  In  this 
case,  the  new  critical  pressures  for  Models  17  ano  20  should  have  been  0.2589  of 
the  critical  pressures  of  Models  8  and  11  respectively,  i.e.,  16.22  lb/in4  and 
16.15  lb/in2  (compare  with  19.893  lb/in2  and  20.277  lb/in2  of  Table  5). 


A  portion  of  Table  5  was  compiled  using  various  formulae  in  the  open 
literature  to  give  an  estimate  of  the  critical  pressure  at  which  buckling  is  to 
occur,  in  theory,  at  least  in  the  elastic  range.  Furthermore,  it  provides 
estimates  to  the  collapse  of  ring  stiffeners  subject  to  live  pressure.  Note 
that  for  Models  4  through  7  (Model  G)  for  which  the  stiffeners  were  not 


3%irstein,  DTMB  Report  997. 

^Bushnell,  D. ,  "Effect  of  Cold  Bending  and  Welding  on  Buckling  of  Ring- 
Stiffened  Cylinders,"  Journal  of  Computers  and  Structures,  Vcl.  12,  1980, 

pp.  291-307. 
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unae roe signed ,  the  critical  pressure  tor  ring  buckling  by  Table  5  is  lb. 410 
lb/in^.  while  the  experimental  collapse  of  the  sh^ll  25  lb/'in".  The  Southwell 
method  ^  yielaed  19.585  lb/in*"  and  the  Von  Mises"*^  23.091  lb/in  ,  respectively 
for  Models  4  through  7.  The  point  to  be  conveyed  is  that  when  the  approximate 
theoretical  formulae  define  critical  collapse  pressures  of  the  same  relative 
magnitude  as  STAGS,  with  stresses  below  the  yield  point,  buckling  probably 
occurs  between  the  stiffeners  in  a  diamond-shape  fashion.  On  the  other  hand  , 
when  approximate  formulae  suggest  critical  buckling  pressures  for  a  ring 
stiffener  well  below  the  finite  element  predicted  value  with  stresses  below 
yield,  the  probable  mode  of  collapse  would  be  overall  instability,  such  as  in 
Models  8  through  11  and  7  through  20  (Table  5). 


^Moussouros,  NSWC  TR  81-325. 
^^Moussouros ,  NSWC  TR  81-325. 
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CHAPTER  5 
SUMMARY 


This  report  compares  finite  element  static  pressure  predictions  to 
experimental  values  reported  in  the  open  literature.  As  expected,  there  are 
discrepancies  between  these  values.  The  most  obvious  one  appears  to  be  residual 
stresses  or  strains. 

Simply  supported  or  clamped  (as  defined  here)  end  conditions  lead  to 
similar  predictions  for  the  buckling  pressures.  Continuity  conditions,  in 
general,  may  give  the  buckling  pressure  for  a  ring  (oval  mode) ,  but  it  may  be 
preferable  not  to  use  them  except  when  the  stiffeners  are  especially  heavy. ^ 

Considering  the  overall  performance  of  STAGS,  it  is  safe  to  say  that,  with 
some  exceptions,  it  has  had  some  success  in  predicting  critical  pressures, 
especially  in  the  elastic  range.  However,  as  the  thickness  of  pressure  hulls 
increases  and  R/h  decreases,  failure  may  occur  by  plastic  yielding  and 
axisymmetric  collapse  at  a  limit  point.  For  these  cases,  predictions  using 
STAGS  are  substantially  less  successful.  It  must  be  further  stressed,  however, 
that  in  all  llkelihoqd,  even  in  the  plastic  range,  the  presence  of  the  rings, 
viewed  as  an  Imperfection,  will  dominate  buckling. 

In  closing  it  must  be  mentioned  that  it  is  hoped  to  overcome  some  of  the 
limitations  of  the  current  analysis  in  the  future  by  exploiting 

(i)  a  different  analysis  option  of  STAGS 

(ii)  better  experimental  data 

(iii)  other  computer  codes  such  as  BOSOR,  ABAQUS  or  any  other  technique 
that  may  become  available  in  the  meantime. 


43Moussouros,  NSWC  TR  81-325. 
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END  3 


u  *  Axial  displacement  (along  global  longitudinal  axis) 
v  *  Tangential  displacement  (along  local  axis) 
w  ■  Radial  displacement  (along  local  radial  axis) 

Ru  *  Rotation  about  longitudinal  axis  x 
Rv  *  Rotation  about  local  tangential  axis  v 
Rw  ■  Rotation  about  local  radial  axis  w 

FIGURE  1.  FULL  CYLINDER  CIRCUMFERENTIALLY  AND  HALF  AXIALLY  DISPLAYING  CIRCULAR 
ENDS  4  AND  3  AND  ASSUMED  LONGITUDINAL  ENDS  2  AND  4  AT  0-0  AND  d-n 
RESPECTIVELY 
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FIGURE  8.  MODE  1  FOR  MODEL  8  (ALUMINUM  MODEL  G  OF  MIDGLEV  AND  JOHNSON) 
SUBJECT  TO  CONTINUITY  CONDITIONS  AND  W  -  0  AT  END  1 
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FIGURE  12.  MODE  1  FOR  MODEL  11  (ALUMINUM  MODEL  F  OF  MIDGLEY  AND  JOHNSON) 
SUBJECT  TO  CONTINUITY  CONDITIONS  AND  W  ■  0  AT  END  1 
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FIGURE  13. 


MODE  1  FOR  MODEL  12  (STEEL  1/5  MODEL  OF  KINRA)  SUBJECT  TO  CONTINUITY 
CONDITIONS  AND  W  *  0  AT  END  1 
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FIGURE  15.  DEFORMED  SHAPE  OF  MODEL  14  (STEEL  MODEL  BR-4)  AT  S00  Lb/in 
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FIGURE 


TABLE  1.  TEST  MODEL  STRUCTURAL  DETAILS 
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TABLE  2.  TEST  MODEL  MATERIAL  DETAILS 


ACTUAL  DIMENSIONS  USED  IN  COMPUTATIONS 
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TABLE  4.  MODELING  DETAILS 


FEATURE / 
MODEL  NO. 

TOTAL 

DEGREES 

OF 

FREEDOM 

NO.  OF 

NODE  RINGS 
BETWEEN 
FRAMES 

DEGREE  OF 
SYMMETRY 

BOUNDARY 

CONDITION 

AT  END  1* 

BOUNDARY 
CONDITION 
AT  END  2* 

1 

1350 

3 

Half  Axially 
Half  Circum¬ 
ferentially 

RV=RW=0 

ufo 
v¥  o 

RV=RW=0 

U=0 

1A 

1950 

3 

Half  Axially 
Half  Circum¬ 
ferentially 

RV=RW=0 

W=0 

u^o 

RV=RW=0 

U=0 

2 

2646 

3 

Half  Axially 
Full  Circum¬ 
ferentially 

RV=RW=0 

VlfO 

U?*0 

RV=RW=0 

U=0 

3 

1014 

7&3 

Half  Axially 
Half  Circum¬ 
ferentially 

RV=RW=0 

UfO 

U*0 

RV=RW=0 

U=0 

4 

702 

1 

Half  Axially 
Half  Circum¬ 
ferentially 

RV=RW=0 

w=o 

u^o 

RV=RW=0 

U=0 

5 

702 

1 

Half  Axially 
Half  Circum¬ 
ferentially 

RV=RW=0 

w=o 

u?o 

RV=RW=0 

U=0 

6 

702 

1 

Half  Axially 
Half  Circum¬ 
ferentially 

RV=RW=0 

W=0 

U*0 

RV=RW=G 

u=o 

7 

702 

1 

Half  Axially 
Half  Circum¬ 
ferentially 

RV=RW=0 

UfO 

U*0 

RV=RW=0 

U=0 

*At  both  ends  cangential  displacement  V=0  at  90°  (.ana  270J  for  full  moael).  For 
full  Model  No.  2  at  8“0  and  8=180°  W=0.  (In  general  W  is  free,  if  W*0  is 
stated  in  bounaary  conditions). 
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FEATURE  / 
MODEL  NO. 

TOTAL 

DEGREES 

OF 

FREEDOM 

NO.  OF 

NODE  RINGS 
BETWEEN 

FRAMES 

i 

DEGREE  OF 

SYMMETRY 

BOUNDARY 

CONDITION 

AT  END  1* 

BOUNDARY 
CONDITION 
AT  END  2* 

8 

1326 

1 

Half  Axially 
Half  Circum¬ 
ferentially 

V=W=RV=0 

U^O 

RV=RW=0 

U=0 

9 

1326 

1 

Half  Axially 
Half  Circum¬ 
ferentially 

RV=RW=0 

UfO 

U?0 

RV=RW=0 

U=0 

10 

1326 

1 

Half  Axially 
Half  Circum¬ 
ferentially 

RV=RW=0 

WfO 

U^O 

RV=RW=0 

U=0 

11 

1326 

1 

Half  Axially 
Half  Circum¬ 
ferentially 

RV=RW=0 

W=0 

UfO 

RV=RW=0 

U=0 

12 

7650 

3 

Half  Axially 
Half  Circum¬ 
ferentially 

RV=RW=0 

W=0 

U/=0 

RV=RW=0 

U=0 

13 

546 

1 

Half  Axially 
Half  Circum¬ 
ferentially 

1 

RV=RW=0 

W=0 

U^O 

RV=RW=0 

U=0 

14 

546 

1 

Half  Axially 
Half  Circum¬ 
ferentially 

RV=RW=0 

W=0 

Uj^O 

RV=RW=0 

U=0 

15 


858 


Half  Axially 
Half  Circum¬ 
ferentially 


RV=RW=0 

W=0 

U^O 


RV=RW=0 

U=0 


*These  approximate  formulae  employ  elastic  values.  To  incorporate  plasticity  effects  we  need  to  use 
tangent  modulus  etc. 


COLLAPSE  PRESSURE 

COLLAPSE  EXPERIMENTAL  COLLAPSE  (PS I)  PREDICTED  UY 

PRESSURE  COLLAPSE  MODE  EXPERIMENTAL  RING  FORMULA, 

DATA/  PREDICTED  PRESSURE  PREDICTED  COLLAPSE  SOUTHWELL,  VON 
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